Correlation effects in the transport through quantum dots 
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We study the charge and heat transport through the correlated quantum dot with a finite value 
of the charging energy U ^ oo . The Kondo resonance appearing at temperatures below Tk is 
responsible for several qualitative changes of the electric and thermal transport. We show that under 
such conditions the semiclassical Mott relation between the thermopower and electric conductivity 
is violated. We also analyze the other transport properties where a finite charging energy U has a 
significant influence. They are considered here both, in the limit of small and for arbitrarily large 
values of the external voltage eV = [Il — I^r and/or temperature difference Tl—Tr. In particular, 
we check validity of the Wiedemann-Franz law and the semiclassical Mott relation. 
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Recent development in nanoelectronics has renewed 
the interests for studying the correlation effects in sys- 
tems where the localized quantum states coexist and in- 
teract with a sea of itinerant electrons. Previously some 
related ideas have been explored in the solid state physics 
when investigating the screening of magnetic impurities 
hybridized with the conduction band electrons 1]. Cor- 
relations lead there to the Kondo effect Q which shows 
up by a logarithmic increase of the electrical resistance 
at temperatures below Tk- Similar mechanism has been 
later found to be responsible for a superconductivity of 
the heavy fermion compounds [3[. 

In the nano- and mesoscopic systems (such as quan- 
tum dots or carbon nanotubes connected to external 
leads) there can also arise the Kondo effect [|[. This 
has been predicted theoretically 0, H| and observed ex- 
perimentally in measurements of the differential conduc- 
tance G(V) = dI ^y ' 0]. Recent experimental studies of 
the thermoelectric power S = — -¥f i/(y)=o [H provided 
an additional evidence in favor of the Kondo effect at low 
temperatures . 

Various transport properties through the correlated 
quantum dots have been studied in quite a detail [H, [(| 0, 
H, [l(| HH, EH 03 • Authors investigated the quantum 
dots coupled to the normal, superconducting and ferro- 
magnetic leads. In several cases there have been observed 
some qualitatively new effects, for instance a splitting of 
the zero bias resonance under magnetic field 0, [l2j , ab- 
sence of the even-odd parity effects [f(| or the out of equi- 
librium Kondo effect [3, |ll|, [lj] . Kondo effect has been 
found in such nanosystems like the single atom [l5[ , sin- 
gle molecule [l6[ and carbon nanotubes 17 1. Moreover, 
it has been also observed in the quantum dots attached 
to ferromagnetic leads [Tc| where, in principle, transport 
is controlled by the spin degree of freedom. 

Quantum dots are often thought to be promising ob- 
jects for the nanotechnology e.g. in various spintronics 
applications [T^, as the gate s of quantum computers, in 
electron spin entanglers [2(| or efficient energy conver- 
sion devices [2lJ. Therefore a detailed understanding of 
the charge and heat transport through mesoscopic sys- 
tems seems to be of primary importance [22|. Besides 
the differential conductance it is important to analyze on 



equal footing also the thermoelectric properties because 
they are very sensitive to dimensionality as well as the 
electronic spectrum near the Fermi level. In the Coulomb 
blockade regime the thermoelectric coefficient S has been 
shown 23] to oscillate as a function of the gate voltage (a 
characteristic saw-tooth shape). At temperatures of the 
order or smaller than Tk there forms a narrow resonance 
slightly above the Fermi energy. Its appearance is re- 
sponsible for increasing the differential conductance (the 
zero bias anomaly) and simultaneously leads to a change 
of sign of the thermopower from positive (when T > Tk) 
to negative values (for T smaller than Tk)- Both these 
effects are useful for detecting the Kondo effect in the 
quantum dots [H, [24[ ■ 

In most of the experimental measurements it has been 
found that the charging energy (the on-dot Coulomb re- 
pulsion) U is large, yet being finite. Usually it varies be- 
tween 2 and 5 in units of the effective coupling T whose 
meaning will be explained in the Section II. 

It is a purpose of our work here to focus on study- 
ing the transport properties for the quantum dot with 
a finite Coulomb energy U. We analyze the linear re- 
sponse with respect to the bias V and temperature dif- 
ference AT applied across the junction and discuss some 
nonlinear effects in the case of arbitrary V and AT val- 
ues. Transport coefficients are determined with use of the 
non-equilibrium Keldysh Green's functions [25[ . We treat 
the correlations using the equation of motion approach 
[2?| whose justification has been recently thoroughly dis- 
cussed in the Ref. (27j . 

Studying the equilibrium and the non-equilibrium as- 
pects of the electric and thermal transport we shall check 
a validity of the Wiedemann-Franz law and the Mott re- 
lation. The paper is organized as follows. In Section II 
we briefly introduce the model and the approximation 
used for calculating the transport coefficients. Results 
obtained within the linear response theory are presented 
in Section III. In the next Section IV we discuss a non- 
equilibrium case corresponding to arbitrarily large values 
of V and AT. We enclose this paper by a summary and 
some overview of the future prospects. 
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I. FORMULATION OF THE PROBLEM 
A. The model 

To account for the correlation effects we focus on the 
simplest situation where QD is coupled to the normal 
leads as described by the non-degenerate single impurity 
Anderson Hamiltonian 

H = ^ ^00^0^ + ^ e d d t d <y + Un d] n di 



k,/3,<7 



(1) 



Operators Ck^o- ( c u.aa) annihilate (create) electrons in the 
left ((3 — L) or the right h.s. (/3 = R) electrodes with the 
corresponding energies £k/3o- = £ ker — pp measured from 
the chemical potentials shifted by the external voltage 
PL — PR = Operators d a , e?+ refer to the localized 
electrons of the dot which is characterized by a single 
energy level td as well as by the Coulomb potential U. 
The last term in (QJ describes hybridization between the 
localized and itinerant electrons. 

The retarded Green's function of the QD can be ex- 
pressed in a general form 



(2) 



It consists of two contributions 5^ (a;) to the selfenergy, 

\v - 1 2 

where the first part T, (uj) = ^ ^J^J, corresponds to 
the noninteracting case U — 0. Twice of its imaginary 
part is often used as a convenient definition of the effec- 
tive coupling ^p(oj) = 2nJ2k | Vk/3 | 2 5(oj — Ck/s) ■ We assume 
a flat function Tp(w) — Tp for \u\ < D so, that the mi- 
correlated QD is characterized by the Lorentzian density 
of states centered around Sd with the halfwidth + Tr. 

The real many-body problem is encountered in @ via 
additional part £/(a>) whenever U ^ 0. In the extreme 
limit U — ► oo one can use the slave boson approach which 
approximately yields the following Green's function 



where fp(tu) = 1/ [l + exp^i^/fceT)] . For sufficiently 
low temperatures the local density of states Pd{^) — 
— ~hnG r d (u) + i0 + ) develops a narrow Kondo peak at the 
Fermi level as illustrated in figure [TJ 

One of the methods for studying the effects of Coulomb 
interactions U uses a truncation scheme for a set of the 
coupled Green's functions (29J. This, so called, equation 
of motion (EOM) technique is approximate. In particu- 
lar, it is known that such procedure underestimates both 
the width and height of the Kondo resonance [28| . There 
are also other techniques like NCA or the large N ex- 
pansion but all of them have only a limited range of ap- 
plicability depending on the values of U, Vkp as well as 
temperature T. In this work we simply use the EOM to 
explore a qualitative behavior of the transport properties 





M G d {u) = [l-(n d )}/ w - £d -£ k>/3 |V k/3 



FIG. 1: Density of states pd{u) of the QD in the equilibrium 
case. For a finite Coulomb energy U the spectrum consists of 
two Lorentians centered around ed and + U. At sufficiently 
low temperatures T <Tk there appears a narrow Kondo res- 
onance near the Fermi level (here it corresponds to u = 0). 



for the system ((T|). Our results can be eventually quan- 
titatively improved by the more sophisticated methods. 

According to the standard derivation based on the 
EOM and the Keldysh formalism [3(| one obtains the 
following retarded Green's function (2) of the quantum 

dot m 



(3) 



where 

0dM = [u-£dY 
Ex(u;) = ^|K k/3 | 

k/3 



^-£k/3 U3 — U— 2e rf + £ k /3 



S 3 (w) 



k/3 



1 



1 



^-£k£< UJ-U-2Ed + £,k/3 



In figure fl| we show the density of states Pdiyi) com- 
puted for the equilibrium case (V = and Tj, = Tr). The 
right h.s. panel illustrates the spectrum for the case of a 
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finite U which is characterized by two Lorentians, one at 
£<i and another one around Ed + U. The upper Coulomb 
satellite is missing (see the left h.s. panel) when U = oo. 
At sufficiently low temperatures T < Tk there appears a 
narrow Kondo resonance due to the singlet state formed 
from the itinerant electrons and the localized electrons of 
the quantum dot. 



B. Definition of the transport coefficients 

By applying some bias V or temperature imbalance 
Tl 7^ Tr the system is driven out of its equilibrium 
[26| . Effectively there are induced the charge and 
heat currents through the QD. Using the nonequlibrium 
Keldysh Green's functions [25( one can derive the fol- 
lowing Landuer-type expressions for the charge I(V, AT) 
and heat currents Iq(V, AT) [ij] 



I = - 



dojT(ui) [f L (oj)-f R (u;)} p d (w), 



(4) 



dwT(u) (cj-eV) [/l(w)-/r(o;)] p d (w) 



(5) 

We introduced here the symmetrized coupling T(to) — 
2rL(id)Tji(id)/{rL(u) + T_r(cj)] and factor 2 comes from 
the contributions of f and J. electrons in the tunneling. 
We notice that both currents (j4|5|) are convoluted with 
the spectral function of the QD hence the behavior illus- 
trated in figure [1] indirectly affects the transport proper- 
ties. 

In figure [2 we show the differential conductance 
G(V) = dl/dV obtained for temperatures higher and 
lower than Tk- We notice a clear increase of the zero volt- 
age conductance (zero bias anomaly) when T < Tk while 
the other side-peaks in the conductance at \eV\ ~ U 
are not sensitive to temperature. In realistic systems 
there should be even more peaks because of the multi- 
level electronic structure containing a mesoscopic number 
of atoms at QD. We thus emphasize that from these fea- 
tures only the zero bias anomaly is driven by the many- 
body Kondo effect at temperatures T < Tk, usually be- 
ing of the order of hundreds mK [7[ . 

For a further analysis we introduce the following trans- 
port coefficients 



1=0 



-(— 
\AT 



(6) 
(7) 



which describe correspondingly the thermoelectric power 
and the heat conductance ([7]). In the limit of small 
perturbations V and AT the heat conductance is ex- 
pected to obey the Wiedemann- Franz law k/TG(0) = 
7r 2 /3e 2 , however we will show that this relation is, in 
general, not obeyed for temperatures T ~ Tk (as has 
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FIG. 2: The differential conductance G(V) versus the applied 
bias V for temperatures above and below Tk and a finite value 
of the Coulomb interaction U = 50F. 



been independently pointed out by several authors, e.g. 
0). 

The thermoelectric power © is a quantity which is 
very sensitive to a particle-hole symmetry in the density 
of states of QD. For small perturbations V, AT one can 
analytically prove that the thermopower is proportional 
to the energy weighted with respect to the Fermi level /x 



(f£ - M) 
eT 



(8) 



in the range of thermal excitations ksT. Appearance of 
the narrow Kondo peak should thus be accompanied by 
a change of sign in S(T) from the positive (T > T K ) to 
negative (T < Tk) values. Indeed, this behavior has been 
recently observed experimentally Q. In the following 
sections we discuss in a more detail some properties of 
the thermopower for a representative set of U and for 
varying perturbations V, AT ranging from the small to 
the arbitrarily large values. 



II. SMALL PERTURBATIONS V AND AT 

For infinitcsimally small perturbations V and AT one 
can expand the currents Q O up to linear terms 



(9) 
(10) 



These coefficients Ly of the linear response theory can 
be determined from the corres pon ding correlation func- 
tions in the Kubo formalism [31] . The zero bias con- 
ductance is then given by G(0) = ^ in, the Seebeck 
coefficient by S — hp k 12 - and the thermal conductiv- 

J el La 

ity by k = yj- ^L 2 2 — ■ Coefficients Ljj can also be 
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FIG. 3: Temperature dependence of the zero bias differen- 
tial conductance G(V = 0) obtained from the linear response 
calculations using the coefficient L\\ (fTTjl . 



derived directly from equations Q O and they become 
functional of the QD density of states via 



L H — T 



L\2 — —r- 

h 



L 



22 



T2 

h 



dtu T(lu) 

>g 
o 

dus T{ui) u 



dm 



d u) 



dT 
dT 



(11) 
(12) 



Pd(Lo) (13) 



where Pd(w) and derivatives of the Fermi function are to 
be computed for V = 0. From the Onsager relation we 
have L 2 i = Li2- 

From the equation (fTTj) we conclude that at very low 
temperatures the zero bias conductance G(0) is propor- 
tional to the density of states pd{oJ — 0). Appearance 
of the Kondo automatically enhances the zero bias con- 
ductance as has been indeed observed experimentally Q • 
Moreover, for T — > it reaches the unitary limit value 
2e 2 /h. With an increase of temperature there activated 
the scattering processes involving some higher energy sec- 
tors [33. 

The linear response relations (|1HI13[) are restricted to 
the external voltage V and/or temperature difference AT 
being small (in terms of T). In the remaining part of 
this section we give a brief summary of our numerical 
computations while the next section will be devoted to 
arbitrary large perturbations. 



Numerical results 



Using the linear response relations (|1HI13[) we explored 
numerically the transport properties for several values of 
T, U, Ed keeping a fixed number of charge on the QD 
(rid) = 0.5. We assumed the flat coupling Tp^ui) = T for 
energies \u\ < D and Tp(w) — elsewhere. 

Figure [3] shows the temperature dependence of the 
zero bias conductance. We notice that G(0) decreases 
monotonously versus temperature. At low temperatures 
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FIG. 4: Temperature dependence of the thermal conductance 
k(T) obtained in the linear response theory. 
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FIG. 5: The Wiedemann Franz relation obtained in the linear 
response theory for the same values of U as in figure [4] 



T < O.ir, conductance decreases because of a gradual 
disappearance of the Kondo peak. For higher tempera- 
tures beyond the Kondo regime (T ~ T) there is a bit 
of a plateau and then again there starts an exponential 
decrease with respect to T. This behavior seems to be 
universal and the Coulomb interaction U has merely a 
marginal effect on it. 

On contrary, the effect of Coulomb interactions does 
show up in the temperature dependence of the thermal 
conductivity. In figure we can see2]that at small temper- 
atures k(T) exponentially increases versus T. For T ~ T 
there appears one maximum (when the values of U is 
small so that the Lorentzians presented in figure[T]overlap 
with one another) and then the other maximum occurs 
at T ~ U (for sufficiently large U). The heat current 
is thus enhanced either when transport occurs through 
the low energy sector (the Lorentzian around Ed) or it 
has an additional contribution by activating the high en- 
ergy sectors (the Lorentzian around Ed + U). We checked 
that the Kondo state has no influence on k(T). The mu- 
tual relation between k(T) and the G(0) is summerized 
in figured] 
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FIG. 6: Temperature dependence of the thermopower S(T) 
obtained in the linear response theory for several values of 
the Coulomb interaction: U = 50T (dotted line), U = 25T 
(dashed line) and U = 17. 5F (solid curve). 



Figure [6] illustrates the Seebeck coefficient S(T). At 
low temperatures T < Tk the thermopower becomes 
negative due to the Kondo resonance slightly above the 
Fermi level (figure[T]) . A finite value of the Coulomb inter- 
action U leads to a partial flattening of the minimum in 
the Kondo regime. At higher temperatures S(T) changes 
sign and considerably increases. For T ~ U there occurs 
an other change of sign (negative S(T)) because of acti- 
vating the high energy sector (the upper Forentian in fig- 
ure[T]). Similar behavior has been independently reported 
for the periodic Anderson model [33] . At extremely large 
temperatures a magnitude of the thermopower asymp- 
totically scales as \S(T) \ oc T^ 1 . 



B. The Mott relation 

Within the linear response limit (i.e. for small pertur- 
bations V and AT) the following Mott formula [3J] 



S M (T) =- — 
K ' 3 



7T 2 k%T <91nG(0) 



Of i 



(14) 



is often used in experimental studies of the thermopower 
S(T). Dependence of the zero bias conductance G(V — 
0) on the chemical potential can be in practice masured 
by varrying the gate voltage Vq. Since the gate voltage 
shifts the energy levels of the QD one can assume that 
- d lnG(0)/d/z ~ d \nG(0)/dV G . 

For the realistic multilevel QDs the Mott relation ([Fi]) 
implies for the thermopower to have a sawtooth shape 
as a function of Vg [231 ] - Similar property can be ob- 
served even in the single level QD as a consequence of 
the Coulomb blockade [35]. The dashed lines in figure 
[7] show the variation of the thermopower obtained from 
the Mott relation (fT4")) for temperatures higher (the upper 
panel) and lower than the Kondo temperature Tk (the 
bottom panel). In the Kondo regime we notice a clear 
discrepancy between the Mott formula and the results 
computed directly from ([J) and ©■ 
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FIG. 7: The thermopower S(T) obtained directly from the 
linear response theory (solid lines) and using the Mott relation 
tUt (dashed lines) . 



In order to clarify a limittcd range of applicability for 
the semi-classical Mott relation we briefly recollect some 
basic constraints used for deriviation of (|T4j) as has been 
also independently stressed by other outhors [35]. The 
Landauer-type expression Q can be expanded with re- 
spect to the small perturbations V and AT. If one ne- 
glects the temperature dependence of pd{^>) (which ob- 
viously is not valid for T ~ Tk) then such expansion is 
carried out only in the Fermi distribution functions fp(u)) 



d u 



uj — yL 
T 



(Tp-T) 



.(15) 



Substituting /i^ = /i+ ^eV and /ir = fx— \eV together 
with (fl"5"|) into the equation ((4]) we search for such a bias 
V which compensates the current induced by the tem- 
perature imbalance AT = T^ — Tr. Using the definition 
([6|) we finally obtain (3f| 



S(T) = 



f-oa - ^)r(w)pd(^) 


f a/Ml 




dtu 





(16) 



For a quantitative determination of (JT6J) one can ap- 
ply the Sommerfeld expansion J_ du> f (u>)Y (u>) ~ 

dwY{u) + ^-(k B T) 2 Y' (n) which is valid only at low 
temperatures. Within the lowest order estimation one 
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FIG. 8: The thermopower S(T) computed within the linear 
response theory in a range of the gate voltage Vg where the 
Kondo effect has the most strong influence. 



FIG. 9: The thermopower S(T) = Q determined beyond 

the linear response theory in the Kondo regime for Ed = — 3r 
and U = 50r. 



finally gets 

S(T) = 



f_ fc|r 91np d (w = 0) k=0 
3 e 



dpi 



(17) 



because limy- 



g/M 



_ — S(u—fi). According to (fTTj) 
the zero bias conductance lim^^o G(V = 0) = const x 
Pd(uJ = 0)\ T=o hence we finally arrive at the relation (fl4|) 



S(T) 



7T 2 jfe|T 91nG(F = 0)| T=o 



(18) 



It is now clear that deviation from the Mott relation 
shown in figure [7] are due to the temperature dependence 
of the QD spectrum Pd{^) when the Kondo peak gradu- 
ally builds in. 

In figure [5] we present the qualitative changes of the 
thermopower caused by the Kondo effect. For high tem- 
peratures T > Tk we see a piece of the saw-tooth be- 
havior [23] . Upon lowering the temperatures to T < Tk 
the Kondo resonance is formed in the density of states 
(see figure [IJ and through the relation (fT6|) this affects a 
sign of the thermopower. Such phwical effects have been 
recently observed experimentally [8| . 

It is interesting to notice that the qualitative changes 
of the thermopower occur on both teeth of the saw. We 
explain in the appendix that this is related to the particle 
or hole type of the Kondo resonance. 



III. BEYOND THE LINEAR RESPONSE 

We also studied the charge (|4]) and heat currents ([5]) go- 
ing beyond the limit of small bias V and temperature dif- 
ference AT. Figure [5] shows the differential conductance 
computed numericaly for Tl = Tr. When temperatures 
Tp differ from one another the differential conductance 
has still qualitatively the similar behaviour with only the 
zero bias anomaly being gradually smeared for increasing 
AT. 
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FIG. 10: Variation of the thermal conductance k(T) as a func- 
tion of the applied temperature difference AT at low temper- 
ature region (upper panel) and at higher temperatures (the 
bottom panel). 



To account for nolincar effects of the other transport 
quantities we determied the thermopower and thermal 
conductivity directly from the definitions ([B][7]). Follow- 
ing [24| , we solved numerically the integral equation 



I(V, AT) = 



(19) 



for the fixed temperature difference AT determining the 
bias V(AT). The thermopower is then given by the ratio 
([5]) . To compute the thermal conductance we determined 
the heat current Iq (V, AT) and substituted it to ([7]) using 
V(AT). 

In figure [HI we show influence of the finite tempera- 
ture difference AT on the thermopower S(T) as function 
of the average temperature T = (Tl + Tr)/2. Our re- 
sults prove that the nonlinear effects play a considerable 
role. In general, by increasing the temperature differ- 
ence AT we notice that the system acquires such prop- 
erties which are characteristic for the high temperature 
behaviour shown previously in figure [6l 

Thermal conductance k is a physical quantity which 
is rather weakly sensitive to any particular structure of 
the low energy excitations. In the linear response limit 
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FIG. 11: The Wiedemann-Franz relation beyond the linear 
response limit determined for several values ol Tl with the 
Coulomb energy U = 25F. On abscissa we put the everage 
temperature defined as T = Tl + 



AT 
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it is known to vanish for T — > as well as for the op- 
posite limit T — > [3, H3|- The maximum value of k 
occurs usually at temperatures T ~ T/fce (see figure H]). 
This situation changes radically if the applied tempera- 
ture difference AT is large. Obviously, k(T) gradually 
increases as a function of AT (see the figure [T0|) but this 
dependence is irregular. We further illustrate it in figure 
[TT] where we plot the Wiedemann-Franz ratio versus the 
average temperature T for several fixed values of Tl- 
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APPENDIX 

In this appendix we briefly explain under what con- 
ditions there can arise the Kondo effect in the system 
described by Hamiltonian (TT]). For this purpose let us 
first explore the particle hole transformation 



d a -> d~l 
dt — > cL 



Ck/3cr 
C+ 



C k/3er 



Ck/3cr 



Substituting these operators into the Hamiltonian (TT]) we 
obtain 



H = £ fa 



k,p,a 



i^c^fia + ^edd+do- + Un d ^n dl 

CkfJadv + Vkff°k/3<rd<r) + COnst 



k,p,a 



where 



IV. SUMMARY 

We studied the influence of correlations on various 
properties of the charge and heat transport through the 
quantum dot connected to the normal leads. In the 
regime of linear response theory (for the small applied 
bias V or temperature difference AT) we find the qualita- 
tive changes caused strictly by appearance of the Kondo 
resonance at low energies. In particular, it leads to the 
zero bias anomaly and to the change of sign of the ther- 
mopower. Moreover, the Coulomb interactions have ad- 
ditional effects by activation of the high energy channels 
in the QD spectrum for the transport of charge and heat. 
Presence of such high energy states manifests by the non- 
monotonous variation of the thermopower (figure [6]) and 
thermal conductance (figured]). 

Appearance of the Kondo effect has a further strong in- 
fluence on deviation from the semiclassical Mott-relation. 
The Mott formula (H]) fails for T ~ Tk because of the 
temperature dependence of the spectral function. 

The Wiedemann-Franz law seems to be rather well pre- 
served unless temperatures are small, for incresing tem- 
peratures we observe a systematic departure from the 
well known Wiedemann-Franz ratio The nonlinear effects 
studied within the Landauer formalism show an addi- 
tional effect on the transport properties for practically 
all transport quantities. This issue has been previously 
pointed also by some other authors [24j |. 



6c/3 

U 



£k/3 

e d -U 



= U 

= -v, 



kp 



(A.l) 
(A.2) 
(A.3) 
(A.4) 



and const = U + 2e d + X)k p a £k/3<r ■ We thus notice 
that the Hamiltonian JTJ) preserves its structure under 
the particle-hole transformation and simultaneously the 
model parameters are scaled as given in (|A.HIA.4j) . This 
property will be useful for us when determining condi- 
tions necessary for the Kondo effect to appear. 

The standard way to specify when electron of the QD 
is perfectly screened by the mobile electrons is obtained 
using the Schrieffer-Wolf transformation. This pertur- 
bative method eliminates the hybridization terms V^p to 
linear terms. In the present case one derives the following 
effective superexchange coupling [3a ] 



J 



u\v k[ 



■P\ 



(e d + U)(-e d ) 



(A.5) 



In the equilibrium case (i.e. for V = 0) this formula (|A.5P 
means that antiferromagnetic interactions arise when the 
energy level e d of the QD is located slightly below the 
chemical potential (we assumed here ^ — 0). This situa- 
tion is schematically illustrated in figure 1. 

The other possibility for the Kondo state to appear 
corresponds to holes rather than particles. In a simple 
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way one can verify that under the particle hole transfor- 
mation the superexchange coupling transforms to 

J= fr|*k.fll a , . (A.6) 
{-i d )(i d + U) 

This result could be expected on more general grounds 



due to the fact that the XY model is invariant on the 
particle hole transformation. We hence conclude that 
the Kondo effect should be present for the small negative 
value of id, in other words for Ed located slightly above 
U. We would like to emphasize that this is not an artifact 
depending on the approximations used in this work. 
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